Hardware quantum gate

ABSTRACT

A hardware quantum gate for running quantum algorithms in a very fast manner exploits the fact that a large number of multiplications required by an entanglement operation of the quantum algorithm provides a null result since only one component per row of the entanglement matrix U F  is not a null. The entanglement operation generates an entanglement vector by permuting pairs of opposite components of a linear superposition vector, depending on the value assumed by the function f. More specifically, if function f is null in correspondence to the vector identified by the first (leftmost) n qubits in common with the two n+1 qubit vectors, in which a pair of opposite components of the superposition vector is referred to, then the corresponding pair of components of the entanglement vector is equal to that of the superposition vector, otherwise it is the opposite.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present application is a continuation-in-part of U.S. patent application Ser. No. 10/615,446 filed Jul. 8, 2003, the entire contents of which are incorporated herein by reference.

FIELD OF THE INVENTION

The present invention relates to quantum algorithms, and more precisely, to a hardware quantum gate that performs, in particular, the entanglement operation of quantum algorithms.

BACKGROUND OF THE INVENTION

Quantum algorithms are global random searching algorithms based on the principles, laws and effects of quantum mechanics. They are used for controlling a process or for processing data in a database, and more specifically, for controlling a process that may include search-of-minima intelligent operations.

In a quantum search, each design variable is represented by a finite linear superposition of initial states, with a sequence of elementary unitary steps manipulating the initial quantum state |i> (for the input) such that a measurement of the final state of the system yields the correct output. Usually, three principle operators, i.e., linear superposition (coherent states), entanglement and interference, are used in the quantum search algorithm.

For a better understanding, a brief description of quantum search algorithms is provided. The problems solved by quantum algorithms may be stated as follows:

Input A function f: {0, 1}^(n) →{0, 1}^(m) Problem Find a certain property of f The structure of a quantum algorithm is outlined by a high level representation in the schematic diagram of FIG. 1.

The input of a quantum algorithm is always a function f from binary strings into binary strings. This function is represented as a map table, which defines for every string its image. Function f is first encoded into a unitary matrix operator U_(F) depending on f properties. This operator calculates f when its input and output strings are encoded into canonical basis vectors of a Complex Hilbert Space: U_(F) maps the vector code of every string into the vector code of its image by f.

$\quad\begin{matrix} {\mspace{146mu}{{BOX}\mspace{14mu} 1\text{:}\mspace{14mu}{UNITARY}\mspace{14mu}{MATRIX}\mspace{14mu} U_{F}}} \\ {{A\mspace{14mu}{squared}\mspace{14mu}{matrix}\mspace{14mu} U_{F}\mspace{14mu}{on}\mspace{14mu}{the}\mspace{14mu}{complex}\mspace{14mu}{field}\mspace{14mu}{is}\mspace{14mu}{unitary}}\mspace{14mu}} \\ {{if}\mspace{14mu}{its}\mspace{14mu}{inverse}\mspace{14mu}{matrix}\mspace{14mu}{coincides}\mspace{14mu}{with}\mspace{14mu}{its}\mspace{14mu}{conjugate}\mspace{14mu}{transpose}\text{:}} \\ {\mspace{256mu}{U_{F}^{- 1} = U_{F}^{=}}} \\ {A\mspace{14mu}{unitary}\mspace{14mu}{matrix}\mspace{14mu}{is}\mspace{14mu}{always}\mspace{14mu}{reversible}\mspace{14mu}{and}\mspace{14mu}{preserves}\mspace{14mu}{the}\mspace{20mu}{norm}} \\ {{of}\mspace{14mu}{{vectors}.}} \end{matrix}$

When the matrix operator U_(F) has been generated, it is embedded into a quantum gate G, a unitary matrix whose structure depends on the form of matrix U_(F) and on the problem to be solved. The quantum gate is the core of a quantum algorithm. In every quantum algorithm, the quantum gate acts on an initial canonical basis vector (the same vector can always be chosen) to generate a complex linear combination (called a superposition) of basis vectors as the output. The superposition contains all the information to answer the initial problem.

After the superposition has been created, a measurement takes place to extract this information. In quantum mechanics, measurement is a non-deterministic operation that produces as output only one of the basis vectors in the superposition. The probability of every basis vector being the output of a measurement depends on its complex coefficient (probability amplitude) in entering a complex linear combination.

The segmental action of the quantum gate and of the measurement forms the quantum block. The quantum block is repeated k times to produce a collection of k basis vectors. In measuring a non-deterministic operation, these basic vectors would not be necessarily identical and each one of them will encode a piece of the information needed to solve the problem. The last part of the algorithm includes interpretation of the collected basis vectors to get the right answer for the initial problem with a certain probability.

The behavior of the encoder block is described in the detailed schematic diagram of FIG. 2. Function f is encoded into matrix U_(F) in three steps.

Step 1: The map table of function f:{0,1}^(n)→{0,1}^(m) is transformed into the map table of the injective function F:{0,1}^(n+m)→{0,1}^(n+m) such that: F(x ₀ , . . . , x _(n−1) , y ₀ , . . . , y _(m−1))=(x ₀ , . . . , x _(n−1) , f(x ₀ , . . . , x _(n−1))⊕(y ₀ , . . . , y _(m−1)))  (1)

$\quad\begin{matrix} {\mspace{146mu}{{{BOX}\mspace{14mu} 2\text{:}\mspace{14mu}{XOR}\mspace{14mu}{OPERATOR}}\mspace{14mu} \oplus}} \\ {{{The}\mspace{14mu}{XOR}\mspace{14mu}{operator}\mspace{14mu}{between}\mspace{14mu}{two}\mspace{14mu}{binary}\mspace{14mu}{strings}\mspace{14mu} p\mspace{14mu}{and}\mspace{14mu} q}\mspace{14mu}} \\ {{of}\mspace{14mu}{length}\mspace{14mu} m\mspace{14mu}{is}\mspace{14mu} a\mspace{14mu}{string}\mspace{14mu} s\mspace{14mu}{of}\mspace{14mu}{length}\mspace{14mu} m\mspace{14mu}{such}\mspace{14mu}{that}\mspace{14mu}{the}\mspace{14mu} i\text{-}{th}\mspace{14mu}{digit}} \\ {{of}\mspace{14mu} s\mspace{14mu}{is}\mspace{14mu}{calculated}\mspace{14mu}{as}\mspace{14mu}{the}\mspace{14mu}{exclusive}\mspace{14mu}{OR}\mspace{14mu}{between}\mspace{14mu}{the}\mspace{14mu} i\text{-}{th}\mspace{14mu}{digits}} \\ {{of}\mspace{14mu} p\mspace{20mu}{and}\mspace{20mu} q\text{:}} \\ {\mspace{205mu}{p = \left( {p_{0},\ldots\mspace{11mu},p_{n - 1}} \right)}} \\ {\mspace{205mu}{q = \left( {q_{0},\ldots\mspace{11mu},q_{n - 1}} \right)}} \\ \left. {s = {{p \oplus q} = \left( {{\left( {p_{0} + q_{0}} \right)\mspace{11mu}{mod}\mspace{14mu} 2},\ldots\mspace{14mu},{\left( {p_{n - 1} + q_{n - 1}} \right)\mspace{11mu}{mod}\mspace{11mu} 2}} \right)}} \right) \end{matrix}$

The need to deal with an injective function comes from the requirement that U_(F) is unitary. A unitary operator is reversible, so it cannot map two different inputs in the same output. Given that U_(F) is the matrix representation of F, F is supposed to be infective. If the matrix representation of function f is directly used, a non-unitary matrix could be obtained since f could be non-injective. Injectivity is thus fulfilled by increasing the number of bits and considering function F instead of function f. Function f can always be calculated from F by putting (y₀, . . . , y_(m−1))=(0, . . . , 0) in the input string and reading the last m values of the output string.

Step 2: Function F map table is transformed into U_(F) map table, following the following constraint: ∀sε{0,1}^(n+m) :U _(F)[τ(s)]=τ[F(s)]  (2) The code map τ:{0,1}^(n+m)→c² ^(n+m) (c² ^(n+m) is the target Complex Hilbert Space) is such that:

$\begin{matrix} \left. {{\left. {{\tau(0)} = {\begin{pmatrix} 1 \\ 0 \end{pmatrix} = \left| 0 \right.}} \right\rangle\mspace{14mu}{\tau(1)}} = {\begin{pmatrix} 0 \\ 1 \end{pmatrix} = \left| 1 \right.}} \right\rangle & (3) \\ \left. {{\tau\left( {x_{0},\ldots\mspace{11mu},x_{n + m - 1}} \right)} = {{{\tau\left( x_{0} \right)} \otimes \ldots \otimes {\tau\left( x_{n + m - 1} \right)}} = \left| {x_{0}\ldots\mspace{11mu}{x_{n + m - 1}}_{\mspace{11mu}}} \right.}} \right\rangle & \; \end{matrix}$

BOX  3:  VECTOR  TENSOR  PRODUCT⊗ The tensor product between two vectors of dimensions h and k is a tensor product of dimension hk, such that: $\left. {\left. {❘x} \right\rangle \otimes} \middle| y \right\rangle = {{\begin{pmatrix} x_{1} \\ \ldots \\ x_{h} \end{pmatrix} \otimes \begin{pmatrix} y_{1} \\ \ldots \\ y_{k} \end{pmatrix}} = \left. \begin{pmatrix} {x_{1}y_{1}} \\ \ldots \\ {x_{1}y_{k}} \\ \ldots \\ {x_{h}y_{1}} \\ \ldots \\ {x_{h}y_{k}} \end{pmatrix}\Rightarrow\begin{matrix} \underset{\_}{{Physical}\mspace{14mu}{interpretation}\text{:}} \\ \text{If a component of a complex vector is interpreted as the probability amplitude of a system being in a given state (indexed by the component number, the tensor product  between two vectors describes the joing probability amplitude of two systems being in a  joint state.} \end{matrix} \right.}$

EXAMPLES Vector Tensor Products

${\left( {0,0} \right)\overset{\mspace{20mu}\tau\mspace{20mu}}{\rightarrow}{\begin{pmatrix} 1 \\ 0 \end{pmatrix} \otimes \begin{pmatrix} 1 \\ 0 \end{pmatrix}}} = {\begin{pmatrix} 1 \\ 0 \\ 0 \\ 0 \end{pmatrix} = \left. 00 \right\rangle}$ ${\left( {0,1} \right)\overset{\mspace{20mu}\tau\mspace{20mu}}{\rightarrow}{\begin{pmatrix} 1 \\ 0 \end{pmatrix} \otimes \begin{pmatrix} 0 \\ 1 \end{pmatrix}}} = {\begin{pmatrix} 0 \\ 1 \\ 0 \\ 0 \end{pmatrix} = \left. 01 \right\rangle}$ ${\left( {1,0} \right)\overset{\mspace{20mu}\tau\mspace{20mu}}{\rightarrow}{\begin{pmatrix} 0 \\ 1 \end{pmatrix} \otimes \begin{pmatrix} 1 \\ 0 \end{pmatrix}}} = {\begin{pmatrix} 0 \\ 0 \\ 1 \\ 0 \end{pmatrix} = \left. 10 \right\rangle}$ ${\left( {1,1} \right)\overset{\mspace{20mu}\tau\mspace{20mu}}{\rightarrow}{\begin{pmatrix} 0 \\ 1 \end{pmatrix} \otimes \begin{pmatrix} 0 \\ 1 \end{pmatrix}}} = {\begin{pmatrix} 0 \\ 0 \\ 0 \\ 1 \end{pmatrix} = \left. 11 \right\rangle}$

Code τ maps bit values into complex vectors of dimension 2 belonging to the canonical basis of C². Besides using a tensor product, τ maps the general state of a binary string of dimension n into a vector of dimension 2 ^(n) for reducing this state to the joint state of the n bits composing the register. Every bit state is transformed into the corresponding 2-dimensional basis vector, and then the string state is mapped into the corresponding 2^(n)-dimensional basis vector by composing all bit-vectors through a tensor product. The tensor product is thus the vector counterpart of a state conjunction. Basis vectors are denoted using the notation |i>. This notation is taken from the description of quantum mechanics.

Step 3: U_(F) map table is transformed into U_(F) using the following transformation rule: [U _(F)]_(ij)=1

U _(F) |j>=|i>  (4) which can be easily understood when considering vectors |i> and |j> as column vectors. Associating these vectors to the canonical basis, U_(F) defines a permutation map of the identity matrix rows. In general, row |j> is mapped into row |i>. This rule will be illustrated in detail in an example quantum algorithm: Grover's algorithm.

The core of the quantum block is the quantum gate, which depends on the properties of matrix U_(F). The scheme in FIG. 3 gives a more detailed description of the quantum block. The matrix operator U_(F) in FIG. 3 is the output of the encoder block represented in FIG. 2. Here, it becomes the input for the quantum block.

This matrix operator is first embedded into a more complex gate: the quantum gate G. Unitary matrix G is applied k times to an initial canonical basis vector |i> of dimension 2 ^(n+m). Every time, the resulting complex superposition G|0 . . . 01 . . . 1> of basis vectors is measured, and one basis vector |x_(i)> is produced as a result. All the measured basis vectors {|x₁>, . . . , |x_(k)>} are collected together. This collection is the output of the quantum block.

The intelligence of such algorithms is in the ability to build a quantum gate that is able to extract the information necessary to find the required property of f and to store it into the output vector collection. The structure of the quantum gate for every quantum algorithm will be discussed in detail, observing that a general description is possible. To represent quantum some special diagrams called quantum circuits are going to be used.

An example of a quantum circuit, relative to the so called Deutsch-Jozsa's quantum algorithm, is illustrated in FIG. 4. Every rectangle is associated to a matrix 2^(n)×2^(n), where n is the number of lines entering and leaving the rectangle. For example, the rectangle marked U_(F) is associated to matrix U_(F). Typically, matrix H represents a Hadamard rotation

$\begin{matrix} {H = {\frac{1}{\sqrt{2}}\begin{bmatrix} 1 & 1 \\ 1 & {- 1} \end{bmatrix}}} & (5) \end{matrix}$

Quantum circuits provide a high-level description of the gate, and using some of the transformation rules listed in FIGS. 5 a-5 f, it is possible to compile them into the corresponding gate-matrix.

$\quad\begin{matrix} {\mspace{76mu}\underset{\_}{\;{{BOX}\mspace{14mu} 4\text{:}\mspace{14mu}{MATRIX}\mspace{14mu}{TENSOR}\mspace{14mu}{{PRODUCT}\mspace{11mu} \otimes}}}} \\ {{{The}\mspace{14mu}{tensor}\mspace{14mu}{product}\mspace{14mu}{between}\mspace{14mu}{two}\mspace{14mu}{matrices}\mspace{14mu} X_{n \times m}\mspace{14mu}{and}{\;\mspace{11mu}}Y_{h \times k}}\;} \\ {{is}\mspace{14mu} a\mspace{14mu}({block})\mspace{14mu}{matrix}\mspace{14mu}\left( {n \cdot h} \right) \times \left( {m \cdot k} \right)\mspace{14mu}{such}\mspace{14mu}{that}\text{:}} \\ {\;{{X \otimes Y} = {{\begin{bmatrix} {x_{11}Y} & \cdots & {x_{1m}Y} \\ \cdots & \cdots & \cdots \\ {x_{n1}Y} & \cdots & {x_{n\; m}Y} \end{bmatrix}\mspace{14mu}{with}\mspace{14mu} X} = \begin{bmatrix} x_{11} & \cdots & x_{1m} \\ \cdots & \cdots & \cdots \\ x_{n1} & \cdots & x_{n\; m} \end{bmatrix}}}} \end{matrix}$

Example Matrix Tensor Product

$\overset{\_}{\underset{\_}{{\begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} \otimes \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix}} = {\begin{bmatrix} {1 \cdot \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix}} & {2 \cdot \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix}} \\ {3 \cdot \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix}} & {4 \cdot \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix}} \end{bmatrix} = \begin{bmatrix} 5 & 6 & 10 & 12 \\ 7 & 8 & 14 & 16 \\ 15 & 18 & 20 & 24 \\ 21 & 24 & 28 & 32 \end{bmatrix}}}}$ It will be clearer how to use these rules when the first example of a quantum algorithm will be discussed in greater detail below.

The decoder block has the function to interpret the basis vectors collected after the iterated execution of the quantum block. Decoding these vectors means to retranslate them into binary strings and interpreting them directly if they already contain the answer to the starting problem. Alternatively, the vectors may be used as coefficient vectors for an equation system to get the searched solution. This part will not be discussed in any greater detail because it is relatively straightforward to understand by those skilled in the art.

Because of the particular importance of the Grover's quantum algorithm in the realization of controllers and data search algorithms in databases, a brief description of the Grover's algorithm will be given below. Grover's problem is stated as follows:

Input A function f: {0, 1}^(n)→{0, 1} such that ∃x∈{0, 1}^(n): (f(x) = 1

∀y∈{0, 1}^(n): x ≠ y

f(y) = 0) Problem Find x

In the Deutsch-Jozsa's algorithm there are two classes of input functions and it must be determined what class the input function belongs to. In this case the problem is almost identical in its form, even if it is more difficult because now we are dealing with 2^(n) classes of input functions (each function of the kind described forms a class).

The diagram of the Grover's algorithm is depicted in FIG. 6, and the gate equation is Φ=[(D _(n) {circle around (×)}I)·U _(F)]^(h)·(^(n+1) H)  (6) Operator D_(n) is called a diffusion matrix of order n and is responsible for interference in this algorithm. The diffusion matrix is defined as follows:

D_(n) |0 . . . 0> |0 . . . 1> . . . |i> . . . |1 . . . 0> |1 . . . 1> |0 . . . 0> −1 + 1/2^(n−1) 1/2^(n−1) . . . 1/2^(n−1) . . . 1/2^(n−1) 1/2^(n−1) |0 . . . 1> 1/2^(n−1) −1 + 1/2^(n−1) . . . 1/2^(n−1) . . . 1/2^(n−1) 1/2^(n−1) . . . . . . . . . . . . . . . . . . . . . . . . |i> 1/2^(n−1) 1/2^(n−1) . . . −1 + 1/2^(n−1) . . . 1/2^(n−1) 1/2^(n−1) . . . . . . . . . . . . . . . . . . . . . . . . |1 . . . 0> 1/2^(n−1) 1/2^(n−1) . . . 1/2^(n−1) . . . −1 + 1/2^(n−1) 1/2^(n−1) |1 . . . 1> 1/2^(n−1) 1/2^(n−1) . . . 1/2^(n−1) . . . 1/2^(n−1) −1 + 1/2^(n−1) Grover's algorithm may be implemented in routines for searching a desired item in a set, by representing in vector form each item of the set forming an input set of vectors, and applying a Grover's algorithm to this set of vectors. The output vector represents the desired item.

The implementation of the Grover's algorithm clearly implies the calculation of several vector products. In fact, all qubits must be multiplied by matrix H, then by the entanglement matrix U_(F) and all qubits but the latter must be multiplied by matrix D_(n).

These multiplications could be carried out via software, but it is quite evident that the number of qubits of a quantum algorithm is very critical in terms of computational speed. In fact, referring to the scheme in FIG. 6, the addition of only one qubit doubles the dimensions of the matrices. Thus, the number of elements (and products) increases exponentially.

A method of performing the superposition operation of the Grover's or Deutsch-Jozsa's quantum algorithm over an input set of vectors is disclosed in European Patent No. 1267304, which is assigned to the current assignee of the present invention. This method exploits the fact that any rotated vector, obtained performing the Hadamard rotation (on an input vector) contemplated by the superposition operation of these quantum algorithms can be easily encoded in a binary vector. Therefore, the successive tensor product of the rotated vectors for generating linear superposition vectors can be carried out by logic gates. This fact allows a noticeable time saving since logic gates are very fast.

However, this is not sufficient to significantly speed up running these quantum algorithms because the entanglement matrix U_(F) is a 2^(n+1)×2^(n+1) square matrix, which implies a considerable computational weight both in the Grover's algorithm as well as in the Deutsch-Jozsa's algorithm.

Differently from other quantum algorithms, in the Grover's algorithm it is possible to iterate h times the entanglement and interference operations until the best solution is reached. An example of evolution of the Grover's algorithm with n=3 is given in FIG. 7 a, in which basis vector and superposition, entanglement and interference output vectors are reported in order. Several iterations of entanglement and interference operations produce a better distribution of probability amplitudes.

Each value is a component on the output vector referred to a basis of vectors of n+1 qubits. There are couples or pairs of values having an opposite sign, referred to vectors of the basis having in common the first (leftmost) n qubits. For example, the values 0.625 and −0.625 are respectively referred to as vectors |0110> and |0111>, respectively. Each pair of elements having an opposite sign represents the probability amplitude of a certain element of the database. For the considered example, the value 0.625 is the probability of the element associated to vector |011> after 3 iterations (h=3).

The algorithm may be iterated as far as a certain quantity is to be minimized, calculated as a function of the components of the output vector, and is smaller than a certain pre-established value. For instance, this quantity can be the Shannon entropy:

$\begin{matrix} {{S(h)} = {- {\sum\limits_{k = 1}^{2^{n + 1}}\;{{{q_{k}(h)}}^{2}\log{{q_{k}(h)}}^{2}}}}} & (7) \end{matrix}$ where q_(k)(h) is the k-th component of the output vector Q taken after h iterations.

The components of the output vector obtained after h=15 iterations are represented in FIG. 7 b. From FIG. 7 b it is clear that the element of the database to be found is associated with vector |011>, and after 15 iterations the Grover's quantum algorithm will find it with a probability of about 0.69.

From the above discussion it is evident that the problem of the large number of computations is even more crucial in the Grover's algorithm than in the Deutsch-Jozsa's algorithm because multiplication by the entanglement matrix U_(F) and the interference matrix D_(n){circle around (×)}I might be repeated many (h) times to output the best result.

SUMMARY OF THE INVENTION

In view of the foregoing background, an object of the present invention is to provide a hardware quantum gate for running quantum algorithms in a very fast manner.

The present invention exploits the fact that a large number of multiplications required by the entanglement operation of quantum algorithms give a null result, because only one component per row of the entanglement matrix U_(F) is not a null. The entanglement operation generates an entanglement vector by permuting the places of pairs of opposite components of a linear superposition vector, depending on the value assumed by the function f(.). More specifically, if function f(.) is null in correspondence to the vector identified by the first (leftmost) n qubits in common with the two n+1 qubit vectors to which a pair of opposite components of which the superposition vector is referred to, then the corresponding pair of components of the entanglement vector is equal to that of the superposition vector. Otherwise, it is the opposite.

Therefore, it is not necessary to calculate the entanglement matrix U_(F) to generate an entanglement vector from a superposition vector, but it is sufficient to copy or invert components of a superposition vector to generate corresponding components of an entanglement vector, depending on the values of the function f(.) processed by the quantum algorithm. This can be easily done using switches input with a pair of components having an opposite value of a superposition vector.

More precisely, the object of the present invention is to provide a quantum gate for running quantum algorithms using a certain binary function defined on a space having a basis of vectors of n qubits, composed of a superposition subsystem carrying out a superposition operation over components of input vectors for generating components of linear superposition vectors referred on a second basis of vectors of n+1 qubits. An entanglement subsystem carries out an entanglement operation over components of the linear superposition vectors for generating components of entanglement vectors. An interference subsystem carries out an interference operation over components of the entanglement vectors for generating components of output vectors.

The entanglement subsystem comprises a command circuit generating a number (2^(n)) of logic command signals encoding the values of the binary function in correspondence to the vectors of the first basis, and an array of multiplexers input with the logic command signals that generate, for each superposition vector, corresponding signals representing components of an entanglement vector.

Each component of the entanglement vector referred to a respective vector of the second basis is equal to the corresponding component of the respective superposition vector if the binary function is null in correspondence to the vector of the first basis formed by the first n qubits of the respective vector of the second basis. Alternatively, each component may be equal to the opposite of the corresponding component of the respective superposition vector if the binary function is not a null in correspondence to the vector of the first basis formed by the first n qubits of the respective vector of the second basis.

BRIEF DESCRIPTION OF THE DRAWINGS

The particular aspects and advantages of the invention will become more evident through the following description of several important embodiments and by referring to the attached drawings, wherein:

FIG. 1 is a block diagram of quantum algorithm in accordance with the prior art;

FIG. 2 is a block diagram of an encoder in accordance with the prior art;

FIG. 3 is a general structure of the quantum algorithm as shown in FIG. 1;

FIG. 4 is a circuit for a Deutsch-Jozsa's quantum gate in accordance with the prior art;

FIG. 5 a shows an example of a tensor product transformation in accordance with the prior art;

FIG. 5 b shows an example of dot product transformation in accordance with the prior art;

FIG. 5 c shows the identity transformation in accordance with the prior art;

FIG. 5 d shows an example of the propagation rule in accordance with the prior art;

FIG. 5 e shows an example of the iteration rule in accordance with the prior art;

FIG. 5 f explains the input/output tensor rule in accordance with the prior art;

FIG. 6 is an example of a circuit forming a Grover's quantum gate in accordance with the prior art;

FIGS. 7 a and 7 b illustrate the evolution of the Grover's quantum algorithm in accordance with the prior art;

FIG. 8 is a detailed view of entanglement and interference subsystems of the quantum gate in accordance with the present invention;

FIG. 9 is a detailed view of section I-a as shown in FIG. 8;

FIG. 10 is an embodiment of the section I-b as shown in FIG. 8;

FIG. 11 is an embodiment of an adder as shown section I-c of FIG. 8;

FIG. 12 is a detailed view of the digital section of the quantum gate in accordance with the present invention; and

FIG. 13 is one possible embodiment of the level shifters as shown in FIG. 12.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

The quantum gate of the invention is suitable for a fast running quantum algorithm applied over a set of input vectors, such as for example, decision making or data search routines based on the Deutsch-Jozsa's algorithm or the Grover's algorithm. It is composed of a superposition subsystem carrying out a linear superposition, an entanglement subsystem carrying out an entanglement operation, and an interference subsystem carrying out an interference operation according to the quantum algorithm to be implemented.

An essential feature of the quantum gate of the invention includes the fact that the entanglement subsystem does not multiply a superposition vector for the entanglement matrix U_(F), but generates components of an entanglement vector simply by copying or inverting respective components of the superposition vector depending on values of the function f(.).

This allows a relevant reduction of the number of multiplications with respect to known methods, and can be carried out using multiplexers. For the sake of simplicity, the hardware quantum gate of the invention will be described in reference to the Grover's quantum algorithm for n=3, though what will be discussed may be easily repeated for other quantum algorithms (in particular, the Deutsch and the Deutsch-Jozsa's algorithm) whose entanglement operation includes a permutation of components of a superposition vector.

A hardware quantum gate of the invention suitable for running the Grover's algorithm with any number of iterations is substantially composed of two parts. Part I: (analog) Calculation step-by-step of the output values. This part is divided in the following subsections: I-a: Entanglement; I-b and I-c: Interference. Part II: (digital) Entropy evaluation, storage of vectors for iterations and output display. This part provides also the first basis of vectors.

The analog part of this scheme for a three-qubits quantum gate is depicted in FIG. 8. A command circuit HB14 generates eight command signals Vc1, . . . , Vc8 each representing a value of the function f(.) on a respective vector of the first basis. An array of multiplexers HB13 is input with voltage signals 011, . . . , 082 representing the sixteen components of a linear superposition vector and generates signals Vo1, . . . , Vo8 representing only the even or the odd components of an entanglement vector. Let us suppose that these signals represent the odd components of an entanglement vector.

Section I-a is present in every hardware quantum gate of the invention, irrespective of the quantum algorithm to be implemented. Each multiplexer, depicted in FIG. 9 is input with a pair of components (011 and 012 for example) that are referred to vectors of the second basis having the first (leftmost) 3 qubits in common (|0000> and |0001>), and generates a respective component (Vo1) which is equal to 011 if Vc1 represents a null value of function f(.) on vector |000>, or to 012 if Vc1 represents a non-null value.

Sections I-b and I-c depicted in FIG. 8 are specifically designed for Grover's quantum algorithms. The presence of tensor products in the interference operation, whose number increases dramatically with the dimensions, forms a critical point.

Sections I-b and I-c of FIG. 8 of the quantum gate allow the interference operation of Grover's quantum algorithm to be quickly carried out. It has been noticed that the matrix D_(n){circle around (×)}I has the following properties: odd columns (or rows, because D_(n){circle around (×)}I is symmetric) have non-zero odd components and even columns have non-zero even components. The value of all non-zero components, except for the i^(th) component of i^(th) column (diagonal elements), is ½^(n−1). The components on the up-left down-right diagonal of the matrix differ from the other non-zero components because they are decreased by 1. The variable G* in an entanglement vector, and the output vector of the quantum algorithm V=(D_(n){circle around (×)}I)G* involves only a suitable weighted sum of components of G*. The value ½^(n−1) depends only from the number n of qubits.

From the above analysis, the generic element v_(i) of V can be written as follows as a function of components g*_(i) of the entanglement vector G*:

$\begin{matrix} {v_{i} = \left\{ \begin{matrix} {{\frac{1}{2^{n - 1}}{\sum\limits_{j = 1}^{2^{n}}\; g_{{2j} - 1}^{*}}} - {g_{i}^{*}\mspace{14mu}{for}\mspace{14mu} i\mspace{14mu}{odd}}} \\ {{\frac{1}{2^{n - 1}}{\sum\limits_{j = 1}^{2^{n}}\; g_{2j}^{*}}} - {g_{i}^{*}\mspace{14mu}{for}\mspace{14mu} i\mspace{14mu}{even}}} \end{matrix} \right.} & (16) \end{matrix}$ Therefore, to calculate a component v_(i) of the output vector it is sufficient to calculate a weighted sum of even

$\left( {\frac{1}{2^{n - 1}}{\sum\limits_{j = 1}^{2^{n}}\; g_{2j}^{*}}} \right)$ or odd

$\left( {\frac{1}{2^{n - 1}}{\sum\limits_{j = 1}^{2^{n}}\; g_{{2j} - 1}^{*}}} \right)$ components of the entanglement vector and to subtract from it the corresponding component g*_(i) of the entanglement vector.

An adder HB25, depicted in greater detail in FIG. 10, sums these components with a certain scale factor that depends only on the number n of qubits (which is 0.25 for n=3), and generates a signal SQ representing the sum of the odd (or even) components of the entanglement vector. Finally, an array of adders I-c generates signals 01, . . . , 08 representing the odd (or even) components of an output vector by subtracting the signals Q1, . . . , Q8 from the scaled sum SQ.

A hardware quantum gate of the invention for carrying out a Deutsch-Jozsa's quantum algorithm does not have the sections I-b and I-c of FIG. 8 for performing the interference operation. In this case, the interference subsystem is substantially composed of an array of adders, each calculating a respective component 01, . . . , 08 of the output vector as a linear combination of components Vo1, . . . , Vo8 of a corresponding entanglement vector.

It is not necessary to calculate all components of the entanglement or output vectors, because the odd components of any vector are always opposite to the even components. For this reason entanglement and interference operations are carried out only on the odd or even components. The other components are calculated by inverting the first ones.

The adders of section I-c may be formed, for instance, as depicted in FIG. 11. A voltage of 2.5V is added to transpose the value inside the range [0÷5]V of A/D converters that interface with Part II (digital part).

To provide better performances, a different range of probability amplitudes is chosen. The presence of converters, whose range is [0÷5]V, suggests adoption of a range of [−2.5÷2.5] instead of [−1/√{square root over (2)}÷1/√,{square root over (2)}]. This implies that the output is not normalized to 1 but to 2*2.5²=12.5.

Because of the necessity of different types of operations (storing values, entropy evaluation and stopping iteration) for performing the Grover's algorithm, a microprocessor unit has been chosen as the core of Part II. For instance, an appropriate microprocessor unit may be the commercial device CPLD XC95288XL of STMicroelectronics, as depicted in FIG. 12.

An analog/digital converter, which may be for example the A/D converter ADC0808 of STMicroelectronics, receives signals representing components 01, . . . , 08 of the output vector and produces a corresponding binary string D1, . . . , D8. The microprocessor unit XC95288XL receives this string and calculates the Shannon entropy.

Summarizing, the microprocessor unit performs the following functions:

-   -   1. Drives correctly the A/D converter.     -   2. Acquires eight digital values and evaluates Shannon entropy         S.     -   3. Compares S with a fixed threshold.     -   4. If S<threshold, stops iteration and sends the results to LED         Matrix display; otherwise sends the results to a digital/analog         converter.     -   5. Provide initial condition of superposed basis vectors.

A display may also be connected to the CPLD to display the result. If Shannon entropy is not a minimum, the binary string has to be re-converted in an analog signal by a digital/analog converter for it to feedback into the entanglement subsystem I-a. In the preferred embodiment of FIG. 12, the digital/analog converter is the commercial device AD7228 of STMicroelectronics.

The level shifters HB3, . . . , HB10, which can be formed as depicted in FIG. 13, re-translate the values of IN1, . . . , IN8 into the range [−2.5÷2.5] and re-obtain all components from odd components. The iterations are carried out very fast, according to the CPLD frequency (20÷30 MHz). This fact suggests that very high performances could be reached even with a greater number of qubits, resulting in a very efficient search in a large database. 

1. A processing device for implementing a quantum gate for running a quantum algorithm for searching elements in a database using a binary function having a first basis of vectors of n qubits, the processing device comprising: a superposition subsystem for performing a superposition operation over components of input vectors for generating components of linear superposition vectors for a second basis of vectors of n+1 qubits; an entanglement subsystem for performing an entanglement operation over components of the linear superposition vectors for generating components of entanglement vectors, said entanglement subsystem comprising a command circuit for generating a plurality of logic command signals encoding values of the binary function corresponding to the first basis of vectors, and an array of multiplexers, each being driven by a respective logic command signal and receiving as input a plurality of signals representing components of a linear superposition vector corresponding to the second basis of vectors having the first n qubits in common, and outputting for each superposition vector corresponding signals representing components of an entanglement vector, each component of the entanglement vector corresponding to a respective vector of the second basis of vectors and being equal to the corresponding component of the respective superposition vector if the binary function is null in correspondence to the vector of the first basis of vectors formed by the first n qubits of the respective vector of the second basis of vectors, or the opposite of the corresponding component of the respective superposition vector if the binary function is not a null in correspondence to the vector of the first basis of vectors formed by the first n qubits of the respective vector of the second basis of vectors; and an interference subsystem for performing an interference operation over components of the entanglement vectors for generating components of output vectors representing the elements searched in the database.
 2. A quantum gate according to claim 1, wherein the quantum algorithm comprises a Grover's quantum algorithm, and wherein said interference subsystem comprises: an adder receiving as input voltage signals representing even or odd components of an entanglement vector, and generating a sum signal representing a weighed sum with a scale factor of the even or odd components of the entanglement vector; and an array of adders each receiving as input a respective signal representing an even or odd component of an entanglement vector and the sum signal, and generating a signal representing an even or odd component of the output vector as a difference between the sum signal and the signal representing an even or odd component of the entanglement vector.
 3. A quantum gate according to claim 2, further comprising an elaboration subsystem comprising: an analog/digital converter receiving as input the signals representing the even or odd components of the output vector and for converting the input signals into a digital string; and a microprocessor unit receiving as input the digital string for performing the following calculating a quantity to be minimized that is associated with the odd or even components of the output vector, comparing the quantity to be minimized with a threshold and stopping the Grover's algorithm if the quantity is smaller than the threshold, otherwise commanding another iteration, and generating an output digital string representing components of the output vector.
 4. A quantum gate according to claim 3, wherein said elaboration subsystem further comprises: a digital/analog converter receiving as input the output digital string, and generating output signals corresponding to the odd or even components of the output vector; and an array of level shifters, each being input with a respective output signal, and generating a pair of voltage signals in a desired voltage range representing opposite components of a new superposition vector to be input to said entanglement subsystem.
 5. A quantum gate according to claim 4, wherein each level shifter comprises: an adder that subtracts a predetermined voltage from a respective output signal for generating an odd component of the new superposition vector; and an inverter receiving as input the odd component for generating the corresponding even component.
 6. A quantum gate according to claim 3, wherein the quantity to be minimized comprises a Shannon entropy.
 7. A quantum gate According to claim 1, wherein the quantum algorithm comprises a Deutsch-Jozsa's quantum algorithm, and wherein said interference subsystem comprises an array of adders each receiving as input signals representing even or odd components of an entanglement vector, and generating a signal representing a corresponding even or odd component of the output vector as a linear combination of the signals representing the even or odd components of the entanglement vector.
 8. A processing device for implementing a quantum gate for running a Grover's quantum algorithm for searching elements in a database using a binary function having a first basis of vectors of n qubits, the processing device comprising: a superposition subsystem for performing a superposition operation over components of input vectors for generating components of linear superposition vectors for a second basis of vectors of n+1 qubits; an entanglement subsystem for performing an entanglement operation over components of the linear superposition vectors for generating components of entanglement vectors, said entanglement subsystem comprising a command circuit for generating a plurality of logic command signals encoding values of the binary function corresponding to the first basis of vectors, and an array of multiplexers, each being driven by a respective logic command signal and receiving as input a plurality of signals representing components of a linear superposition vector corresponding to the second basis of vectors having the first n qubits in common, and outputting for each superposition vector corresponding signals representing components of an entanglement vector; and an interference subsystem for performing an interference operation over components of the entanglement vectors for generating components of output vectors representing the elements searched in the database.
 9. A quantum gate according to claim 8, wherein each component of the entanglement vector corresponding to a respective vector of the second basis of vectors is equal to the corresponding component of the respective superposition vector if the binary function is null in correspondence to the vector of the first basis of vectors formed by the first n qubits of the respective vector of the second basis of vectors, or the opposite of the corresponding component of the respective superposition vector if the binary function is not a null in correspondence to the vector of the first basis of vectors formed by the first n qubits of the respective vector of the second basis of vectors.
 10. A quantum gate according to claim 8, wherein said interference subsystem comprises: an adder receiving as input voltage signals representing even or odd components of an entanglement vector, and generating a sum signal representing a weighed sum with a scale factor of the even or odd components of the entanglement vector; and an array of adders each receiving as input a respective signal representing an even or odd component of an entanglement vector and the sum signal, and generating a signal representing an even or odd component of the output vector as a difference between the sum signal and the signal representing an even or odd component of the entanglement vector.
 11. A quantum gate according to claim 10, further comprising an elaboration subsystem comprising: an analog/digital converter receiving as input the signals representing the even or odd components of the output vector and for converting the input signals into a digital string; and a microprocessor unit receiving as input the digital string for performing the following calculating a quantity to be minimized that is associated with the odd or even components of the output vector, comparing the quantity to be minimized with a threshold and stopping the Grover's algorithm if the quantity is smaller than the threshold, otherwise commanding another iteration, and generating an output digital string representing components of the output vector.
 12. A quantum gate according to claim 11, wherein said elaboration subsystem further comprises: a digital/analog converter receiving as input the output digital string, and generating output signals corresponding to the odd or even components of the output vector; and an array of level shifters, each being input with a respective output signal, and generating a pair of voltage signals in a desired voltage range representing opposite components of a new superposition vector to be input to said entanglement subsystem.
 13. A quantum gate according to claim 12, wherein each level shifter comprises: an adder that subtracts a predetermined voltage from a respective output signal for generating an odd component of the new superposition vector; and an inverter receiving as input the odd component for generating the corresponding even component.
 14. A quantum gate according to claim 11, wherein the quantity to be minimized comprises a Shannon entropy.
 15. A processing device for implementing a quantum gate for running a Deutsch-Jozsa's quantum algorithm for searching elements in a database using a binary function having a first basis of vectors of n qubits, the processing device comprising: a superposition subsystem for performing a superposition operation over components of input vectors for generating components of linear superposition vectors for a second basis of vectors of n+1 qubits; an entanglement subsystem for performing an entanglement operation over components of the linear superposition vectors for generating components of entanglement vectors, said entanglement subsystem comprising a command circuit for generating a plurality of logic command signals encoding values of the binary function corresponding to the first basis of vectors, and an array of multiplexers, each being driven by a respective logic command signal and receiving as input a plurality of signals representing components of a linear superposition vector corresponding to the second basis of vectors having the first n qubits in common, and outputting for each superposition vector corresponding signals representing components of an entanglement vector; and an interference subsystem for performing an interference operation over components of the entanglement vectors for generating components of output vectors representing the elements searched in the database.
 16. A quantum gate according to claim 15, wherein each component of the entanglement vector corresponding to a respective vector of the second basis of vectors is equal to the corresponding component of the respective superposition vector if the binary function is null in correspondence to the vector of the first basis of vectors formed by the first n qubits of the respective vector of the second basis of vectors, or the opposite of the corresponding component of the respective superposition vector if the binary function is not a null in correspondence to the vector of the first basis of vectors formed by the first n qubits of the respective vector of the second basis of vectors.
 17. A method for operating a processing device comprising a superposition subsystem, an entanglement subsystem and an interference subsystem for running a quantum algorithm for searching elements in a database using a binary function having a first basis of vectors of n qubits, the method comprising: performing a superposition operation with the superposition subsystem over components of input vectors for generating components of linear superposition vectors for a second basis of vectors of n+1 qubits; performing an entanglement operation with the entanglement subsystem over components of the linear superposition vectors for generating components of entanglement vectors, the entanglement operation comprising generating a plurality of logic command signals encoding values of the binary function corresponding to the first basis of vectors, and driving an array of multiplexers by respective logic command signals and applying as input to the array of multiplexers a plurality of signals representing components of a linear superposition vector corresponding to the second basis of vectors having the first n qubits in common, and outputting for each superposition vector corresponding signals representing components of an entanglement vector, each component of the entanglement vector corresponding to a respective vector of the second basis of vectors and being equal to the corresponding component of the respective superposition vector if the binary function is null in correspondence to the vector of the first basis of vectors formed by the first n qubits of the respective vector of the second basis of vectors, or the opposite of the corresponding component of the respective superposition vector if the binary function is not a null in correspondence to the vector of the first basis of vectors formed by the first n qubits of the respective vector of the second basis of vectors; and performing an interference operation with the interference subsystem over components of the entanglement vectors for generating components of output vectors corresponding to the elements searched in the database.
 18. A method according to claim 17, wherein the quantum algorithm comprises a Grover's quantum algorithm; and wherein the interference subsystem comprises: an adder receiving as input voltage signals representing even or odd components of an entanglement vector, and generating a sum signal representing a weighed sum with a scale factor of the even or odd components of the entanglement vector; and an array of adders each receiving as input a respective signal representing an even or odd component of an entanglement vector and the sum signal, and generating a signal representing an even or odd component of the output vector as a difference between the sum signal and the signal representing an even or odd component of the entanglement vector.
 19. A method according to claim 18, further comprising: receiving as input the signals representing the even or odd components of the output vector, and converting the input signals into a digital string; and performing the following based upon the digital string calculating a quantity to be minimized that is associated with the odd or even components of the output vector, comparing the quantity to be minimized with a threshold and stopping the Grover's algorithm if the quantity is smaller than the threshold, otherwise commanding another iteration, and generating an output digital string representing components of the output vector.
 20. A method according to claim 19, further comprising: generating output signals corresponding to the odd or even components of the output vector; and using an array of level shifters, each being input with a respective output signal, and generating a pair of voltage signals in a desired voltage range representing opposite components of a new superposition vector to be input to an entanglement subsystem performing the entanglement operations.
 21. A method according to claim 20, wherein each level shifter comprises: an adder that subtracts a predetermined voltage from a respective output signal for generating an odd component of the new superposition vector; and an inverter receiving as input the odd component for generating the corresponding even component.
 22. A method according to claim 19, wherein the quantity to be minimized comprises a Shannon entropy.
 23. A method according to claim 17, wherein the quantum algorithm comprises a Deutsch-Jozsa's quantum algorithm; and wherein the interference subsystem comprises an array of adders each receiving as input signals representing even or odd components of an entanglement vector, and generating a signal representing a corresponding even or odd component of the output vector as a linear combination of the signals representing the even or odd components of the entanglement vector. 